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A detection or nondetection of primordial non-Gaussianity by using the cosmic mi- 
crowave background radiation (CMB) offers a way of discriminating inflationary scenarios 
and testing alternative models of the early universe. This has motivated the consider- 
able effort that has recently gone into the study of theoretical features of primordial 
non-Gaussianity and its detection in CMB data. Among such attempts to detect non- 
Gaussianity, there is a procedure that is based upon two indicators constructed from the 
skewness and kurtosis of large-angle patches of CMB maps, which have been proposed 
and used to study deviation from Gaussianity in the WMAP data (see Refs. [I] and 12)1. 
Simulated CMB maps equipped with realistic primordial non-Gaussianity are essential 
tools to test the viability of non-Gaussian indicators in practice, and also to understand 
the effect of systematics, foregrounds and other contaminants. In this work we extend 
and complement the results Refs. [T] and |2] by performing an analysis of non-Gaussianity 
of the high-angular resolution simulated CMB temperature maps endowed with non- 
Gaussianity of the local type, for which the level of non-Gaussianity is characterized by 
the dimensionless parameter /JJ'l'^' ■ 
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1. Introduction 

The statistical properties of the temperature anisotropies of cosmic microwave back- 
ground (CMB) radiation offer a powerful probe of the physics of the primordial 
universe. In particular the study of non-Gaussianity of CMB is a powerful approach 
to probe the origin and evolution of structures in the universe (see, for example, the 
Refs. '3'-^ and references therein). Given the far reaching consequences of a convinc- 
ing detection (or non-detection) of primordial non-Gaussianity for our description 
of the physics of the early universe, it is important to employ different statistical 
tools to quantify its amount, type and the angular scale in order to have information 
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that may be helpful for identifying its causes. Apart from revealing features of non- 
Gaussianity, different statistical estimators can be sensitive to different systematics. 
On the other hand, since one does not expect that a single statistical estimator can 
be sensitive to all possible forms of non-Gaussianity that may be present in CMB 
data, it is important to study the possible deviations from Gaussianity by using 
different statistical tools to identify any non-Gaussian signals in the CMB data. 

Recent analyses of CMB data made with different statistical tools have provided 
indications of either consistency or deviation from Gaussianity (see, e.g., Ref. ?). In 
a recent papeiH we have proposed two new large-angle non-Gaussianity indicators, 
based on skewness and kurtosis of large-angle patches of CMB maps, which pro- 
vide measures of the departure from Gaussianity on large angular scales. We have 
used these indicators to carry out analyses of large-angle deviation from Gaussian- 
ity in both band and foreground-reduced WMAP CMB maps with and without a 
KQ75 mask.flEI Wg found that while non-Gaussianity of the Q, V, and W masked 
maps are consistent with Gaussianity, there is a strong indication of deviation from 
Gaussianity in the K and Ka masked maps. We have also sho wn th at the full-sky 
five-year foreground-reduced internal linear combination as well as the 

harmonic ILC (HILC) and the needlet ILC (NILC)20 maps present a significant 
deviation from Gaussianity.^ 

Simulated CMB maps endowed with assigned primordial non-Gaussianity are es- 
sential tools to test the power and sensitivity of non-Gaussian indicators. The first 
reported simulations of CMB temperature maps with primordial non-Gaussianity 
introduced through a non-Gaussian parameter /nl were given by the WMAP 
team.l^ Subsequently, Liguori et al. produced generalized algorithm that improves 
the computational speed and accuracy,!^ and includes polarization.!^ A set of 
300 temperature and polarization maps with non-Gaus sian ities of the local type 
at the WMAP angular resolution were then produced.l^ More recently, Eisner 
and WandeltP^I presented new algorithm and generated 1 000 high-angular resolu- 
tion simulated non-Gaussian CMB temperature and polarization maps with non- 
Gaussianities of the local type, for which the level of non-Gaussianity is defined by 
the dimensionless parameter /^l^'- 

In this paper, we extend and complement the investigations of Refs|T| and [2] 
(see also the related Ref. [25]) . by using their skewness and kurtosis indicators to 
carry out analyses of Gaussianity of high-angular resolution simulated non-Gaussian 
CMB temperature maps, equipped with non-Gaussianity of local type with different 
amplitude parameters and generated according to the procedure given by 

Eisner and Wandelt^^. 

The structure of the paper is as follows. In Sec. [2] we introduce our non- 
Gaussianity indicators. Section [3] contains the results of applying our statistical 
indicators to the non-Gaussian simulated maps and our main conclusions. 
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2. Indicators and Maps of Non-Gaussianity 

The steps of a constructive way of defining our non-Gaussianity indicators S and 
K, and the associated maps (discrete function s de fined on S"^) from input CMB 
(simulated or real data) maps are the followingPJ2l 

i. Take a discrete finite set of points {j = 1, . . . , A'c} homogeneously distributed on 
the CMB celestial sphere as the centers of spherical caps of a given aperture 
7; and calculate for each cap j the skewness and kurtosis given, respectively, by 

^^^W1^T.(^^-T^)' i^^.^ ^(T, -r,)^-3, (1) 

where A'p is the number of pixels in the j cap, Ti is the temperature at the i ^'^ 
pixel, Tj is the CMB mean temperature of the j cap, and a is the standard 
deviation. Clearly, the whole set of numbers Sj and Kj for each j, obtained 
through this discrete scanning of the CMB sphere, can be viewed as a measure 
of non-Gaussianity in the direction of the center of the cap centered at {9j,(j)j) 
with aperture 7. 

ii. Patching together the Sj and Kj values for all spherical cap j, give discrete 
functions S — S{6, </>) and K — K{d, (j)) defined over the celestial sphere, which 
can be used to measure the deviation from Gaussianity as a function of the 
angular coordinates (6*, The Mollweide projection of skewness and kurtosis 
functions S = S{9, (p) and K — K{9, (p) arc nothing but skewness and kurtosis 
maps, hereafter referred to as 5— map and iiT— map, respectively. 

Clearly, the functions S = S{6, (p) and K = K{d, (p) are functions defined on 
and can be expanded into their spherical harmonics to have their power spectra Si 
and Kg. Thus, for example, for the kurtosis indicator K — K{9, (j)) one has 

00 e 

K{e, 4>)^J2Y1 ^^™(^' ' (2) 

t=0 m=-l 

and can calculate the corresponding angular power spectrum 



21 



which can be used to quantify the amplitude (level) and angular scale of the devi- 
ation from Gaussianity. The power spectrum can also be used to calculate the sta- 
tistical significance of such deviation by comparison with the corresponding power 
spectrum calculated from input Gaussian maps (/nl*^' — 0)- Obviously, similar ex- 
pressions and analyses of the statistical significance can be made for the skewness 
S^S{e,<jy). 

In the next section we shall use the statistical indicators S — S{9, ((>) and 
K — K{9, 4>) to make analyses of non-Gaussianity of the high-angular resolution 
simulated CMB temperature maps endowed with non-Gaussianities of the local 
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type, for which the level of non-Gaussianity is characterized by the dimensionless 

f loca 

/nl 



parameter f^'^'^^ 



3. Main Results and Conclusions 

From the previous section it is clear that in order to calculate the skewness and 
kurtosis functions S = S{6, cj)) and K = K{d, cj)) and the associated S and K 
maps, one ought to have an input CMB map. The input maps used in our anal- 
yses are high-angular resolution simulated CMB temperature maps endowed with 
non-Gaussianities of the local type defined by different values of the dimensionless 
amplitude parameter /^l'''. 

A simulated map with a desired level of non-Gaussianity /nl''' is such that the 
spherical harmonic coefficients are given b; 



a 



L (4) 



where a and a are the linear and non-linear spherical harmonic coefficients 
of the simulated CMB non-Gaussian maps generated in Ref. [24] and are available 
for download^ 
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Fig. 1. Simulated CMB temperature fluctuations maps. The left panel shows a Gaussian map 
(/nl'^' = 0) while the right panel depicts a non-Gaussian map generated with f^{^^^ = 5 000. We 
have taken this value for /nl just as an example that makes the non-Gaussian effects visible by 
naked eye through the comparison between the Gaussian and non-Gaussian maps. Temperatures 
are in mK. 



gives a Gaussian map /^l^' — , while the right panel shows a non-Gaussian maps 
with f^l""^ = 5 000 . The HEALPix resolution parameter for map of Fig. [1] as well 



Figure [T] shows examples of such simulated input CMB maps. The left panel 

floca 
'NL 

.1 c nnn 't^i 

/nl 

as in all the high resolution simulated input CMB maps used in this paper, is 
A^sidc = 512, which corresponds to 3 145 728 pixels. A maximum multipole moment 
^max = 1 024 was also taken. 

In the following we report the results of our analyses made by using 1 000 simu- 
lated as input CMB temperature maps endowed with non-Gaussianities of the local 
type with amplitude parameter /^Iff^' = 0, 100, 1 000, 3 000. 

In our calculations of skewness and kurtosis indicator maps (S— map and 
iC— map), to minimize the statistical noise we have scanned the celestial sphere 

* |http : //planck . mpa-garchlng ■ mpg . de/cmb/f nl-simulat ions [ 
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with spherical caps of aperture 7 = 90°, centered at Nc = 3 072 points on the 
two-sphere homogeneously generated by using HEALPix package!^ 




Fig. 2. Skewness (left) and kurtosis (right) maps calculated from input simulated maps with 
^Nf' ~ ^ 000. Clearly, according to Eq[T]the values in both S and K are dimensionless numbers. 

As an ihustration of typical skewness and kurtosis maps, Fig. [5] shows the Moll- 
weide projection of 5* (left) and K (right) maps generated from an input CMB 
simulated map for /^l^' — 1 000. These maps show spots with higher and lower 
values of S{0, (jy) and K{9, 0), which suggest large-angle dominant components (low 
t) in these maps. We have also calculated similar maps from the simulated input 
maps endowed with non-Gaussianity of local type for the other values of /^l^' that 
we are concerned with in this paper. However, since these maps provide only qual- 
itative information, to avoid repetition we only depict the pair of maps of Figs. [5] 
merely for illustrative purpose. 

In order to obtain quantitative large-angle-scale information of the non- 
Gaussianity of 4 000 calculated S and K maps (obtained from 4 x 1 000 simulated 
input CMB maps generated for = 0, 100, 1 000, 3 000), we have calculated the 

low ^ (£ = 1, • • • , 10) averaged power spectra and Ki, obtained by averaging 
over 1 000 power spectra of S and K maps, calculated for each value of /nl**'- The 
statistical significance of these power spectra is estimated by comparing the val- 
ues of S( and Kg obtained from input maps generated for /^l^' = 100, 1000,3 000 
with the values of the corresponding power spectra and Kg obtained from the 
Gaussian (/nl**' = 0) input simulated map. 

Figure [3] shows the average power spectra of the skewness Sg (left panel) and 
kurtosis Kg (right panel), for £ = 1, • • • , 10 , calculated from input simulated Gaus- 
sian (/nl^' = 0) maps, and from CMB maps equipped with non-Gaussianity of the 
local type for which /jlf^^' = 100, 1000, 3 000. The 95% confidence level, obtained 
from the S and K maps calculated from the Gaussian CMB simulated maps, is 
indicated in this figure by the dashed lineal 

To the extent that the average Se and Kg obtained from input simulated CMB 
maps endowed with = 100 are within 95% Monte-Carlo (MC) average values 

of Se and Kg for /^l''' = 0, Fig. |3] shows that our indicators are not suitable to 

'^For more details on the calculation of S and K maps and associated power spectra we refer the 
readers to Refs. |T]and|2] 
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Fig. 3. Low £ average power spectra of skewness Si (left) and kurtosis (right) Ki calculated from 
the Gaussian (/Jj°£''' = 0) and non-Gaussian {f^l^^^ ^ 0) input simulated CMB maps. The 95% 
confidence level relative to the Gaussian maps is indicated by the dashed line. 



detect this small level of primordial non-Gaussianity of local type in CMB maps. 
However, this figure also shows that they can be effectively employed to detect 
higher level of non-Gaussianity of local type. These results square with the refined 
numerical analysis we shall report in the remainder of this paper. 

To have an overall assessment power spectra Si and Ki, calculated from the 
input simulated non-Gaussian maps equipped primordial non-Gaussianity of local 
type, we have performed a test to find out the goodness of fit for Si and Ke 
multipole values as compared to the expected multipole values obtained from S 
and K maps calculated from Monte-Carlo (MC) statistically Gaussian (/^l**' — 0) 
simulated CMB maps. In each case, this gives a number that quantifies collectively 
the deviation from Gaussianity. For the power spectra Si and Ki we found the values 
given in Table[T]for the ratio x^/dof (dof stands for degrees of freedom) for the power 
spectra calculated from non- Gaussianity of local type with f^l""^ = 100, 1 000, 3 000. 



Table 1. t^st goodness of fit for Si and Ki calcu- 
lated from the maps with different level of non-Gaus- 
sianity as compared with the corresponding expected 
values obtained from MC simulated CMB input maps 
with f^l^' = 0. 



Level of non-Gaussianity 


for Si 


X' for Ki 


f^l'^i = 100 
/Jj-j^'^i = 1000 
f^l^i = 3000 


4.00 X 10"^ 
3.02 X 10 
1.52 X 10^ 


3.10 X 10-2 
5.50 X 10^ 
1.98 X 10"' 



Clearly, the greater are the values for x^/dof the smaller are the x^ probabilities, 
that is the probability that the values of power spectra Si and Ki and the expected 
values of the power spectra of the Gaussian maps agree. Thus, from Table [1] one 
concludes that the maps endowed with /^l**' — 100 present very small level of 
primordial non-Gaussianity, as detected by our indicators, while the non-Gaussinity 
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of the maps with /^l**' — 1 000 and 3 000 are a few orders of magnitude higher and 
suitably detected by both indicators S and K, in particular by the kurtosis indicator 
K. 
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